Propagation of light polarization in a birefringent medium: Exact analytic models 
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Driving the analogy between the coherent excitation of a two-state quantum system and the 
torque equation of motion, we present exact analytic solutions to different models for manipulation of 
polarization in birefringent medium. These models include the one-dimensional model, the Landau- 
Zener model, and the Demkov-Kunike model. We also give an example for robust, broadband 
manipulation of polarization by suitably tailoring the birefringence vector. 
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I. INTRODUCTION 

Various problems in many branches of physics are de- 
scribed by torque equations. These include the textbook 
examples of the classical motion of a charged particle 
in a magnetic field the motion of a point mass un- 
der the Coriolis force as exemplified by the behavior of 
a gyroscope acted on by gravity and the dynamics 
of a spin in a magnetic field [3, |j] . More recent appli- 
cations include the dynamics of optical waveguides [5|, 
polychromatic beam splitters @, sum frequency conver- 
sion techniques in nonlinear optics 0, Q j transformation 
of light polarization in a birefringent medium I9l -ll5l , and 
other nonlinear processes of three wave mixing [l6[ . 

In this paper we explore the analogy between the 
torque equation of motion and the coherent dynamics of 
a two-state quantum system (spin i or a two-state atom) 
to describe analytically the evolution of light polarization 
in a birefringent medium by using the approach of Bloch 
p7| and Fcynman et al. [l^. We use this approach to 
find the propagator for the Stokes polarization vector, 
given the propagator for the respective two-state system. 
This result allows us to find the conditions for some use- 
ful manipulations of the polarization. Then, we derive 
three exact analytic solutions for (i) the one-dimensional 
(resonance) model, (ii) the Landau-Zcner model [l9l.[20j. 
and (ui) the Demkov-Kunike model P.l2l|-[23t. 
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This paper is organized as follows. In section |TT] we 
show the analogy between the torque equation for the 
Stokes polarization vector and coherent excitation of a 
two-state quantum system. In section IIIII we present 
several exact solutions for polarization evolution, in 
which the birefringence vector elements correspond to the 
one-dimensional (resonance) model, the Landau-Zcner 
model, and the Demkov-Kunike model. Within the latter 
two models, we discuss a potentially important example 
for broadband manipulation of polarization. 



II. EVOLUTION OF POLARIZATION 

A. Stokes polarization vector 

Consider a plane electromagnetic wave traveling in 
the z direction through an anisotropic dielectric medium 
without polarization dependent loss (PDL). The light po- 
larization is described by the Stokes polarization vector 
S(z) = [Si{z), 82(2), 83(2)]'^ , which is conveniently de- 
picted on the Poincare sphere shown in Fig. [1] The 
three basic polarizations are described as follows; 

right circular; 8(2) = [0, 0, l]"^; (la) 

left circular; S(z) = [0, 0, -1]^; (lb) 

linear; S(z) = [cos2(/), sin20, 0]"^; (Ic) 

the elliptical polarization is described by the points be- 
tween the poles and the equatorial plane. 

The polarization evolution obeys the following torque 
equation for the Stokes polarization vector jol-flll UBl : 



dz 



S(z) = n{z) X S(z), 



(2) 



FIG. 1: Depiction of the Stokes vector on the Poincare sphere. 



where Jl(z) = [^}i{z),V,2{z),fl3{z)]'^ is the birefringence 
vector of the medium; the direction of fl{z) is given 
by the slow eigenpolarization and its length |r2(2:)| cor- 
responds to the rotary power. The torque equation 
(PI) comprises three coupled linear differential equations 
for the (real) components of the Stokes vector Sk{z) 
{k = 1,2,3). To this end, we use the formal analogy of 
the Stokes polarization equation ([2|) with the Bloch equa- 
tion in quantum mechanics written in Feynman-Vernon- 
Hellwarth's torque form [l8|, and the formal equivalence 
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of the latter to the time-dependent Schrodinger equation 
for a two-state quantum system [23] . Thus the system of 
equations ([2]) can be cast into the form of two coupled 
linear differential equations, 

2i^si(z) = n3{z)si{z) + [ni{z) - m^iz)] s^iz), (3a) 
2i^S2(^) = [ni{z) + in2{z)] si{z) - ^^{z)s2{z), (3b) 

for new complex-valued variables s„(z) (n ~ 1,2). The 
Stokes vector components are bilinear combinations of 
these new variables, 

Si{z) = s\{z)s2{z) + si{z)sl{z), (4a) 

S2{Z) = -I [SI{Z)S2{Z) - Sl{z)s*2{z)] , (4b) 

S,{z) = \s,{zt~\s2{z)\\ (4c) 
Equations ([3]) can be written as a single equation, 

2iAs(z) = H(z)s(z), (5) 

for the vector s(z) — [si(z), S2{z)\^ . Here the Hermitean 
matrix 

H(z) = r!i(z)o-i(z) -f n2{z)a2{z) + n:i{z)(T^{z), (6) 

corresponds to the Hamiltonian in a quantum two-state 
system with (t„ being the Pauli spin matrices. 



where U(zy, z;) is found by using Eqs. dJ) and 0, 

U(z/,Zj) = 

"a2-52-c2 + rf2 -2{ab + cd) 2{ac - bd) 

2{ab-cd) a^-b^ + c^-d^ 2{ad + bc) 
-2{ac + bd) -2{ad - be) + b^ - - d^ 

(9) 

Equation Q allows us to immediately write down the 
conditions for several representative important transfor- 
mations of light polarization. 

(i) Transformation from right circular, S(zi) = 
[0, 0, 1]"^, to left circular polarization, S(z/) = [0, 0, — 1]"^ 
(or vice versa), requires: 

a = 6 = (10) 

which corresponds to a transition probability p = 1 in 
the two-state quantum system. 

(ii) Transformation from right circular, S(zi) = 
[0,0,1]"^, to linear polarization (or vice versa), S(z/) — 
[cos 20, sin 20, 0]-'", requires: 

a2 + 52 = c2 + d2 ^ i (11) 

which corresponds to a transition probability p = ^ in 
the two-state quantum system. 

III. MODELS WITH EXACT ANALYTIC 
SOLUTIONS 



Solution of the Stokes polarization equation in 
terms of a two-state solution 



Because there are a number of analytically exactly sol- 
uble two-state models, we can use the formal correspon- 
dence between the Stokes polarization equation ([2]) and 
the two-state equations ([3]) and derive exact analytic so- 
lutions for the polarization evolution in a birefringent 
medium. Because the "Hamiltonian" ^ is Hermitean, 
the evolution matrix (the propagator) for the variables 
Si{z) and S2{z) is unitary, which is conveniently param- 
eterized by four real parameters a, b, c, d as 



si{zf) 

S2{zf) 



a + lb c + id 
— c + id a — lb 



Sl{Zi) 
S2{Zi) 



(7) 



with a'^ + b^ + c^ + d'^ = 1 . The "transition probability" in 
the two-state problem is p ~ + d'^, whereas + b^ = 
I — p gives the "survival probability" . We express the 
evolution of the Stokes polarization vector in terms of 
the parameters of the two-state propagator. 



Siizf) 

S2{zf) 

S3{zf) 



U(z/,z,) 



Si{zt) 

S2{Zt) 

Ssizi) 



(8) 



Equation ([9]) allows us to use the available exact an- 
alytic solutions for two-state quantum systems to derive 
exact analytic solutions for the evolution of the polariza- 
tion of light propagating through a birefringent medium. 
We select here three such models: the one-dimensional 
(resonance) model, the popular Landau-Zener model 
[H [13 (in its finite version), and the fiexible Demkov- 
Kunike model [22|, [2^ (also in its finite version) . 



A. One-dimensional model 

We assume that the birefringent medium stretches 
from Zi to Zf. The elements of the birefringence vector 
for this model read 



ni{z)^ f{z/L), n2{z) = n^iz) = 0, 



(12) 



where f{z/L) is any function of z for Zi < z < Zf and 
L is a length scale parameter. For example, f{z/L) can 
be chosen in such a way that a = 5 = 0, so we can 
have transformation from left circular to right circular 
polarization, which is known as a half- wave plate [25| . 

This one-dimensional model corresponds to exact res- 
onance in a two-state system, which causes Rabi oscilla- 
tions [26j . and for which the propagator parameters are 
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Distance z (in units L) Distance z (in units l/p ) 



FIG. 2: Transformation of right circular polarization (jlap 
with a birefringence vector Q.{z) = [/(z),0,0]^ for the one- 
dimensional model, where f{z/L) — 1/L. The polarization 
evolution corresponds to Rabi flopping in atomic physics (Eq. 



FIG. 3: Transformation of right circular polarization pap 
with a birefringence vector for the Landau-Zener model with 
flo = P and Zi = -10//3 (Eq. ([TB]), Eq. (O can also be 
applied for z — >■ +oo). 



well known, 



a = cos(A/2), 
d= -sin(A/2), 
6 = c = 0, 



(13a) 
(13b) 
(13c) 



where A{zf,Zi) = J^J f{z / L) dz. Hence the evolution 
matrix for the Stokes vector is 



U(z/,. 



1 

cos(A) -sin(A) 
sin(A) cos(A) 



(14) 



The transformation from right circular to left circular 
polarization is achieved when A = (2fc+ l)7r, /c = 0, 1, 2.., 
i. e. then a = & = 0. An example with the simplest 
function f{z/L) = 1/L, i. e. a constant, is shown in Fig. 
[2j In this case, A = (z/ — Zi)/L and the transformation 
is achieved for — 2;^ = (2fc + 1)ttL. 



For this model, we find the propagator parameters in a 
similar fashion as in |27|. 



F(l - ia^) 

a + 1& = ^ [Aa2(C/)£'-l+ia2(-Cj) 



c + id 



2tt 

_F(l-*a2) 



e'-[Aa^(-C/)Aa2(C.) 



a\/2'K 

- Aa^(C/)Aa^(-G)] 



(16a) 



(16b) 



where D^{z) is the parabolic cylinder (Weber) function 
mill and Ciz) = fize-'T, i. e. 0,/ = /^^z./e'^^; 
a = r2o/2/3. We find the propagator for the evolution 
of the Stokes polarization vector from Eq. ([T5)) by using 
Eq. (HI). 

In the limiting case when L = z/ = — ZiS>l//3 ("sym- 
metric crossing" ) we use the "weak coupling" asymptotic 
of D^(C) (for fixed v and |C| ^ 00) [Ij: 



B. Finite Landau-Zener model 

In this exactly soluble model the birefringence vector 
components are given by 



f)i(z)==r!o, ^^2^=0, l^aW-/?'^, 



(15) 



for Zi < z < Zf, and ftk = (fc = 1,2,3) otherwise. 
These elements correspond to the finite Landau-Zener 
model in atomic physics for the dynamics of the new 
variables si{z) and S2{z), as defined in Eqs. ^ pol [27|. 



DAC) ~ re-^'/4[l + 0(|Cr^)], (largCI < 37r/4), 



(largCI ^3^4), 



(17a) 



(17b) 



and wc find the following approximation of the elements 
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FIG. 4: Transformation of linear polarization (|lcp . = 7r/4, 
with a birefringence vector for the Landau-Zener model with 
0,0 = /? and = (Eq. (fT6l) . Eq. (fT9)) can also be applied 
for z — >■ +oo). 



2 4 6 8 10 
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FIG. 5: Transformation of linear polarization pc|l . (p — 7r/4, 
with a birefringence vector for the Demkov-Kunike model 
with fit, = Ao = l/L, Bo = and Zi = -lOL (Eq. 

(|2ip . Eq. (|23p can also be applied for z — >■ +oo). 



of the propagator [27|: 



c + id= -\f\ 



(18a) 
(18b) 

+ o? ln{l3^L'^) + argr(l - ia^). 

(18c) 



This case is illustrated in Fig. [31 where the exact solution 
based on ([TB]) approaches the asymptotic one ([T7| for 
L = z/ = -Zj > l/,3. 

Another important Hmiting case is when Zi = and 
Zf — L ^ 1/ /3 ("half crossing"). Then, we use D^{0) = 
2'^/20F/r(i^) to simplify Eq. ([TSD [13 



zi (1 2 J 

-'^^^e-^"^/^--cosh(™V2), 



r(i- 2 

jXl_-^a^ 
aV2r(i^) 

y2r(l-ia2; 



(19a) 
(19b) 



ar( 



l — ior 



(19c) 

i^_,„V4+.x-3»^/4 sinh(W/2). 

(19d) 



with X = /32lV4 + a2ln(/3LV2)]. This case is illus- 
trated in Fig. m where the exact solution based on 
approaches the asymptotic one ([T7)) for = and 
2^ = i> 1//3. 



C. Finite Demkov-Kunike model 

We shall derive an analytic solution for birefringence 
vector components given by: 

r2i(z) = 17osech(z/L), (20a) 
172(z)=0, (20b) 
173(2) = Ao + So tanh(z/L), (20c) 

for Zi < z < Zf, and fifc = (fc = 1,2,3), otherwise. 
These elements correspond to the finite Demkov-Kunike 
model for the dynamics of the variables Si{z) and S2{z), 
as defined in Eq. (U) (2^. The solution is derived in a 
similar fashion to the infinite Demkov-Kunike model in 
atomic physics [22|, with the additional assumption of 
finite initial and final conditions. 

The propagator parameters for this model for any Zi 
and Zf are: 



lb 



id = 



-[F,{zf)F*{z,) + U 



22-4S(/3) 



(21a) 



xF,izf)F*iz,)azfy-'az^n 

X \F2{zi)F^{z,)i{zff-^^ -F,{zi)F-,{z,)i{zi)'-^\ 

(21b) 



where a 



/32 



i/3, ^ = -V' 



_BoL/2, 5 = 

_ j^^ J, 



AoL/2, A ^ 

4 + z(<5-/3), 



$0 = \l^^:!,{z')dz\ and $(z,Zi) = f^_VL-i{z')dz' . Fur- 
thermore, Fi(z) = F(A,^, 1^,^(2)) and F2(z) = F(A + 
1 — i^, /i 1 — I', 2 — I/, ^(z)) are hypergeomctric functions, 
where ^(z) = [tanh(z/L) + l]/2 [g^. 
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FIG. 6: Transformation of right circular polarization (|lap 
with a birefringence vector for the Demkov-KunDte model 
with Qo = 1/L, Ao = 1/L, Bo = 1/L, and = (Eq. 
(pTjl . Eq. (f22|) can also be applied for 2 — > +00). 



The formulas for the elements of the propagator for 
any Zi are a generalization of the ones derived in for 
the special case when Zf = — and /3 = (finite Rosen- 
Zener model). Then, the propagator for the evolution 
of the Stokes polarization vector is found from here by 
using Eq. dH). 

When Zi = and Zf — > +00, the elements of the prop- 
agator simplify to [3C| : 



-16 = 
id = 



e'^—Fi 



2(z/ - A - fi) 



(22a) 
(22b) 



where Fi = F{-X, -/i, 1^ - X - fi, F2 = F{1 - A, 1 - 
/i,l + zy — A — /i,i). An example for the manipulation of 
polarization by a birefringence vector with the character- 
istics of the Demkov-Kunike model and Zi = is given 
Fig. E 

The formulas for the elements of the propagator can 
be presented with elementary functions in the case of 
the original Demkov-Kunike model when Zi — —00 and 
Zf +00. Then, the general formula for the elements of 
the propagator can be expressed by [22| : 



16 ^^^S{zf~z,)/L+f3{zf+z,)/L) 



(23a) 
(23b) 



id: 



T{iy - X)T{iy - fi) 

1 — 

r(2 - i')r(iy 



^ - m) 

r(i-A)r(i-AO 



(23d) 



An example for such manipulation of polarization is 
shown in Fig. [S] 



op g 3 

£ u £ 



1.0 





/— X",(z) 








10 -i 


-6 -4 -2 2 4 


6 


8 10 



s 

o 
U 




-6 -4 -2 2 4 6 8 10 
Distance z (in units L) 

FIG. 7: Transformation of right circular polarization (|la[) to 
left circular polarization ()lb|) with a birefringence vector for 
the Demkov-Kunike model with parameters, which satisfy the 
adiabatic condition Qo = 5/1/, Ao = 0, Bo = 5/L (Eq. 



D. Adiabatic limit 

In the limit of large rotary power (which corresponds 
to large pulse area in quantum physics), 

l-Zf 

Global condition: / |r2(z)|dz > 1, (24a) 

J Zi 

Landau-Zener model: fio ^ /? (24b) 
Dcmkov-Kunikc model: VLq > Bq > y ■ (24c) 

the evolution of the Stokes polarization vector becomes 
adiabatic. Additionally, both the Landau-Zener and 
Demkov-Kunike models also involve a "level crossing" 
and they are therefore examples of exactly solvable mod- 
els of rapid adiabatic passage (RAP) l3ll - [33j . 

RAP traditionally provides a well-studied robust and 
efficient method for producing complete population 
transfer between two bound states of a quantum system 
H, [2^ l3ll - [33j . The population change is achieved by 
sweeping the carrier frequency of a laser pulse through 
resonance with a two-state system (parameterized by 
detuning), while simultaneously pulsing the strength of 
the interaction (parameterized as time- varying Rabi fre- 
quency). The sweep of detuning corresponds to a cross- 
ing of the so-called diabatic energy curves. When the 
interaction changes sufficiently slowly, so that the time 
evolution is adiabatic, there will occur a complete popu- 
lation transfer from the ground state to the excited state. 

In terms of manipulation of polarization, RAP is a 
broadband technique for robust adiabatic conversion of 
light polarization and it is therefore applicable to a wide 
range of frequencies. It is also robust to variations in the 
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FIG. 8: Contourplot of S-j,{zj = lOL) where the birefringence 
vector corresponds to the Demkov-Kunike model dynamics 
of the corresponding two-state system with S{zi = — lOL) = 
[0,0,1]^ with parameters: Bo = 5/L. Left circular polariza- 
tion is achieved when S-i{zf) = —1 (Eq. (|21|l '). 



propagation length, the rotary power and other parame- 
ters variations. The robustness of the RAP polarization 
conversion with respect to deviation in the propagation 
length for the Demkov-Kunike model is demonstrated in 
Fig. [71 This figure shows that the Stokes vector is quite 
stable with regard to changes in z for z — > +oo in con- 
trast to the previous figures where it oscillates with z. 
The robustness to variations in the parameters of the 
birefringence vector arc demonstrated in Fig. [8l which 
is a contourplot of S3{zf) with respect to the parame- 
ters Aq and fJp, while the other paremeters are fixed. As 
we can see, right circular polarization is transformed to 
a left circular one for quite large deviations of the two 



parameters. 



IV. CONCLUSION 

In this paper wc have used the analogy between the 
torque equation of motion and the coherent atomic ex- 
citation of a two state system to describe the evolu- 
tion of the Stokes polarization vector in a birefringent 
medium without polarization dependent loss. We have 
given exact analytic solutions for this evolution when 
the birefringent vector corresponds to the dynamics of 
quantum two-state systems, which are described by the 
one-dimensional (resonance) model, the Landau-Zener 
model, and the Demkov-Kunike model. Finally, we have 
demonstrated how this approach could be applied for 
broadband conversion of light polarization in analogy to 
rapid adiabatic passage (RAP) in two-state quantum sys- 
tems. 

This approach can be useful in many other branches 
of physics where we need exact analytic solutions to 
a torque equation. These include the description of 
the Newton equation of motion, the classical motion of 
a charged particle in oscillating magnetic and electric 
fields, e. g. for charged particle confinement, optical 
waveguides, polychromatic beam splitter, sum frequency 
conversion techniques in nonlinear optics, the dynamics 
of a spin in a magnetic field, the behaviour of a gyroscope 
acted on by gravity, etc. 
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